Threshold and Complexity Results for the 
Cover Pebbling Game 



Anant P. Godbole 

East Tennessee State University 

Nathaniel G. Watson 

Washington University in St. Louis 

Carl R. Yerger 

Georgia Institute of Technology 



Abstract 

Given a configuration of pebbles on the vertices of a graph, a pebbling move is 
defined by removing two pebbles from some vertex and placing one pebble on an 
adjacent vertex. The cover pebbling number of a graph, ^{G), is the smallest number 
of pebbles such that through a sequence of pebbling moves, a pebble can eventually 
be placed on every vertex simultaneously, no matter how the pebbles are initially 
distributed. We determine Bose Einstein and Maxwell Boltzmann cover pebbling 
thresholds for the complete graph. Also, we show that the cover pebbling decision 
problem is NP-complete. 
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1 Games 

There are several popular games that involve the movement of "objects" 
along a graph-like structure. These include Mancala, where beads are moved 
along a bent path, and Peg Solitaire, where pegs are moved across a tri- 
angular grid. In each case, some objects are removed from the game after 
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a move is made. Both games proceed according to a set of rules that may 
be seen at http : //www. centralconnector . com/GAMES/mancala.html and 
http://www.mazeworks.com/peggy/ respectively. Mancala is a well-defined 
game between two players, while the solo player in Peg Solitaire pits her- 
self against "nature". This paper concerns the pebbling and cover pebbling 
games, which take place between a highly intelligent Player 1 and a rather 
non- competitive opponent with a limited strategy. There is a version of peb- 
bling, called "pegging" , which is far closer to Peg Solitaire on graphs than are 
our games; see [13]. See also [5] for a chessboard game related to pebbling. 

The focus of our paper, and indeed of all previous research on the subject, 
is on deriving conditions under which Player 1 wins the game, or wins with 
probability that is asymptotic to one, or wins with probability that approaches 
zero as the size of the problem grows to infinity. Let us start with some baseline 
definitions and previously derived facts. 



2 Preliminaries 

Given a connected graph G, distribute t pebbles on its vertices in some con- 
figuration. Specifically, a configuration of weight t on a graph G is a function 
C from the vertex set V{G) to N U {0} such that J2vev(G) C{v) = t. Clearly 
C represents an arrangement of pebbles on V{G). If the pebbles are indis- 
tinguishable, there are ("^^^^) = ("^l^O configurations of t pebbles on n 
vertices. Using quantum mechanical terminology as in [11], we shall call this 
situation Bose Einstein pebbling and posit that the underlying probability 
distribution is uniform, i.e. that each of the ("^1^^) distributions are equally 
likely - should the pebbles be thrown randomly onto the vertices. This is the 
model studied in [8]. Now there is no reason to assume, a priori, that the 
pebbles are indistinguishable. Accordingly, if the pebbles are distinct, we shall 
refer to our process as Maxwell Boltzmann pebbling, in which a random distri- 
bution of pebbles leads to n* equiprobable configurations. Maxwell Boltzmann 
pebbling does not appear to have been studied more than peripherally in the 
literature. 

A pebbling move is defined as the removal of two pebbles from some vertex 
and the placement of one of these on an adjacent vertex. Given an initial 
configuration, a vertex v is called reachable if it is possible to place a pebble 
on it in finitely many pebbling moves. The graph G is said to be pebbleable 
(this is not standard nomenclature) if any of its vertices can be thus reached. 
Define the pebbling number 7r(G) to be the minimum number of pebbles that 
are sufficient to pebble the graph regardless of the initial configuration. The 
pebbhng game may thus be described as follows: Player 2 specifies a distribu- 
tion G and a target vertex v. Player 1 wins the game iff she is able to reach 
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vertex v using a sequence of pebbling moves. The pebbling number of G is the 
smallest number Iq of pebbles so that Player 1 wins no matter what strategy 
Player 2 employs. 

The origin of pebbling is rather interesting and somewhat unexpected. La- 
garias and Saks [14] were considering a way to produce an alternative proof 
to a conjecture of Erdos and Lemke, which Lemke and Kleitman proved in 
1989 [18]. It is known that for any set = {ni, n2, . . . , n^} of q natural num- 
bers, there is a nonempty index set / C {1, . . . , g} such that q\ ^i^jUi. The 
Erdos-Lemke conjecture states that the additional conclusion Yl,iei i^i <lcm 
(g, ni, ■0,2, . . . , Uq) could also be reached ([14]). Unfortunately, Lemke and Kleit- 
man's argument was detailed and contained a considerable amount of case 
analysis. This provoked Lagarias and Saks to invent graph pebbling as a way 
to produce a cleaner proof, since such a proof would follow easily if the peb- 
bling number of the cartesian product of paths was determined. This was 
accomplished in a landmark paper by Chung [4]. One generalization of peb- 
bling, called p-pebbling, was utilized in Chung's proof, and is defined as the 
removal of p pebbles from some vertex and the placement of one pebble on an 
adjacent vertex. It turns out that a greedy-like condition, called the numerical 
pebbling operation for pebbling paths, can be used to prove Chung's theorem 
[14]. In fact, one of the lemmas of Chung's proof actually uses the fact that 
for any set N = {ni,n2, ■ ■ ■ ,nq} of q natural numbers, there is a nonempty 
index set / C {1, . . . , g} such that q\ J2iei ^i- 

SPECIAL CASES: The pebbling number 7r(P„) of the path is 2"-^ Chung [4] 
proved that 7r(Q'^) = 2<^ and 7r(P;") = 2^"-^)"*, where Q'^ is the d-dimensional 
binary cube. An easy pigeonhole principle argument yields 7r(ir„) = n. The 
pebbling number of trees has been determined see [14]. 

One of the key conjectures in pebbling, now proved in several special cases, 
is due to Graham; its resolution would clearly generalize Chung's result for 
m-dimensional grids: 

GRAHAM'S CONJECTURE. The pebbling number of the cartesian product 
of two graphs is no more that the product of the pebbling numbers of the two 
graphs, i.e. 

7r{GnH) < 7r{G)7r{H). 

Structural characteristics of graphs have also been employed to determine the 
pebbling number of specific classes of graphs. For instance, a graph is said 
to be Class if 7r(G') = \G\. Cubes are of Class 0, as are complete graphs, 
but what other families fall in this important class of graphs for which tt is 
as low as it can possibly be? Here are two answers: For graphs of diameter 
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2, if G is 3-connected, that is, the removal of 2 or fewer vertices does not 
disconnect the graph, then G is Class [6]. In fact, if we consider G{n,p), the 
class of random graphs on n vertices where the probability of each particular 
edge being present is a fixed constant p G (0, 1), then almost all such graphs 
are Class [6]. Generalizations of this result to the case where p = Pn ^ 
as n ^ oo are also available; see [10]. Other authors, e.g. [3], have obtained 
general pebbling bounds, while Bukh [2] has proved almost-tight asymptotic 
bounds on the pebbling number of diameter three graphs. 

Another aspect of pebbling that has been explored is the random structure one 
obtains when placing pebbles randomly on graphs. Specifically, we seek the 
probability that a graph G is pebbleable when t pebbles are placed randomly 
on it according to the Bose Einstein or Maxwell Boltzmann scheme. Numerous 
threshold results have been determined in [8] for Bose Einstein pebbling of 
families of graphs such as K^, the complete graph on n vertices; C„, the cycle 
on n vertices; stars; wheels; etc. A threshold result is a theorem of the following 
kind: 

t = tn^ ttn => P(G' = Gn is pcbbleablc) — 1 {n ^ oo) 

t = tn^nbn =^ ^{G = Gn is pcbblcable) — {n ^ oo), 

where we write, for non-negative sequences c„ and c„ ^ dn (or dninCn) if 
Cn/dn — i> oo as n — > oo. Of course, we have reason to feel particularly gratified 
if we can show that a„ = 6„ in a result of this genre. For the families of complete 
graphs, wheels and stars, for example, we know [8] that a„ = 6„ = y/n. In 
many cases, however, the analysis is quite delicate; see [26] and [9] for some 
of the issues involved in finding the pebbling threshold for a family as basic 
as Pn, the path on n vertices. The fundamental reference [1] contains general 
results on the existence of sharp pebbling thresholds for families of graphs. 

A detailed survey of graph pebbling has been presented by Hurlbert [14], and 
it would probably not be an oversimplification to state that most results avail- 
able to date fall in four broad categories: finding pebbling numbers for classes 
of graphs; addressing the issue of when a family of graphs is of class 0; pin- 
pointing graph pebbling thresholds; and seeking to understand the complexity 
issues in graph pebbling [15]. A survey of open problems in graph pebbling 
may be found on Glenn Hurlbert 's website; see http: //math. la. asu. edu/, 
~hurlbert/HurlPebb . ppt. 

The above mini-survey on pebbling notwithstanding, we focus in this paper 
on a variant of pebbling called cover pebbling, first discussed by Crull et al 
[7]. For reasons that will become obvious, we focus only on analogs of the last 
two of the four general directions mentioned above. 

The cover pebbling number \ {G) is defined as the minimum number of pebbles 
required such that it is possible, given any initial configuration of at least A(G) 
pebbles on G, to make a series of pebbling moves that simultaneously reaches 
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eac/i vertex of G. A configuration is said to be cover solvable if it is possible to 
place a pebble on every vertex of G starting with that configuration. Various 
results on cover pebbling have been determined. For instance, we now know 
[7] that \{Kn) = 2n - 1; A(P„) = 2" - 1; and that for trees T„, 

A(r„)= rna^ ^ 2'^^'^-^^\ (1) 

Likewise, it was shown in [16] that A(Q'^) = S'^ and in [24] that A(X^i,...,^^) = 
4ri + 2r2 + . . . + — 3, where ri > r2 > . . . > r„j. The above examples 
reveal that for these special classes of graphs at any rate, the cover pebbling 
number equals the "stacking number", or, put another way, the worst possible 
distribution of pebbles consists of placing all the pebbles on a single vertex. 
The intuition built by computing the value of the cover pebbling number 
for the families Kn, Pn, and T„ in [7] led to Open Question No. 10 in [7], 
which was christened the Stacking Conjecture by students at the Summer 2004 
East Tennessee State University REU. In an exciting summer development, 
participants Annafies Vuong and Ian Wyckoff [22] were able to prove the 

STACKING THEOREM: For any connected graph G, 

X(G) = max V 2'*^^*("'^), 

thereby proving that (1) holds for all graphs. (The stacking theorem was inde- 
pendently proved soon after by Sjostrand [20].) In fact, the key result in Vuong 
and Wyckoff 's paper [22] is really a sufficient condition for a distribution to 
be cover solvable, so further investigations in the theory of (cover) pebbling 
might soon veer, we speculate, in a fifth general direction, namely a study 
of which distributions are (cover) solvable and which are not. Indeed, such a 
research thrust would be most consistent with our description of pebbling as 
a game, and is addressed in Section 5 of this paper. 



3 Mctxwell Boltzmann Cover Pebbling Threshold for 

It is evident that n is the smallest number of pebbles that might suffice to 
cover pebble - in the unlikely event that they happen to be distributed one 

apiece on the vertices. On the other hand, we know that 2n — l pebbles always 
suffice, since X{Kn) = 2n — 1. We seek a sharp cover pebbling threshold that 
is somewhere in between these two extremes, when distinguishable pebbles 
are thrown onto the n vertices of according to the Maxwell Boltzmann 
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scheme. To explain why there is a dramatic increase in the probabiUty of the 

pebbleabihty of Kn at t = (1.5238 )n, we first prove an important auxihary 

result that gives necessary and sufficient conditions for a configuration to be 
cover solvable. Such results are not easy to come by, as we will further see in 
Section 5. 



3. 1 Necessary and Sufficient Conditions for Cover Solvability of 

Let X = Xn^t be the number of vertices on which an odd number of pebbles 
are placed. We will often refer to X as the number of odd stacks. 

Lemma 1 A configuration oft pebbles on the n vertices of is cover solvable 
if and only if 

X + t>2n. (2) 

Proof 1 Let C be any cover solvable configuration. This imphes that after a 
sequence of pebbling moves, each previously uncovered vertex has two pebbles 
associated with it - one on it, and one that was removed from the game. Like- 
wise, any vertex that previously had a non-zero and even number of pebbles 
on it must have at least two pebbles still on it, and any vertex that previously 
had an odd number of pebbles on it must now have at least one pebble on it. 
Thus we must have t > 2E + X (where E represents the number of vertices on 
which an even number (including 0) of pebbles were initially placed), which 
is a condition that is easily seen to be equivalent to (2). 

Conversely, if C is not cover solvable, then we must, after a series of pebbhng 
moves, reach a point where there are still some uncovered vertices; where each 
initially covered vertex in the set E now has exactly two pebbles on it; and each 
vertex in the set X has exactly one pebble on it. It follows that t < 2E + X.D 

Proof 2 Given a configuration C, let Yi,0 < i < t, be the number of vertices 
with i pebbles. Now a vertex can cover exactly i others if and only if it has 
either 2i -|- 1 or 2i -|- 2 pebbles on it. It follows that C is cover solvable iff 

Y^t{Y2i+, + Y,i+2)>Yo, 

i>0 

or, iff 

^(2^ + l)Y2i+i + ^(2^ + 2)^2^+2 > + E ^2i+i + 2 E y2i+2, 

i>0 i>0 i>0 i>0 

i.e.,iSt>2E + X. □ 
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Armed with Lemma 1, we now provide the heuristic reason why we beheve 
there is a sharp cover pcbbUng threshold at t = (1.5238....)n. Given a random 
variable X with expected value IE(X), we will say that X is sharply concen- 
trated around E(X) if X ~ E(X) with high probability (w.h.p.) Assuming 
therefore that X ~ E(X) w.h.p., it makes sense to speculate that Kn is cover 
solvable with high probability whenever E(X) > 2n — t. But X — Y^^=ilj, 
where Ij = 1 (rcsp. 0) if there is an odd (resp. even) stack on vertex j, so that 
linearity of expectation yields 



E(X)=nP(/i = 1) 



j odd 

2 I V n 



"^-fl-^VV (3) 



Thus E{X) >2n-tiS 



and, parametrizing by setting t — An, we see that (4) holds iff 

Since (1 — 2/n)" ~ we see from (5) that a reasonable guess for a threshold 
value of t is Aqu where Aq is the solution of 

^-^exp{-2>l} = | 

or ^40=1.5238 We next make this hunch more precise. 



3.2 Main Result 



Various tools arc used to establish concentration of measure results. Some 
of the more sophisticated techniques employed are the martingale method, 
a.k.a. Azuma's inequality (proved independently and a few years earlier by 
W. Hoeffding), and Talagrand's isoperimetric inequalities in product spaces; 
see [21] for an exposition of both. In this section, however, we establish our 
main result by using a baseline technique, viz. Tchebychev's inequality, also 
known in probabilistic combinatorics circles as the "second moment method." 
Theorem 2 below thus yields the "correct" result, but with a sub-optimal rate 
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of convergence due to the crudeness of the method employed. Tchebychev's 
inequahty states that for any random variable X and for any A > 0, 

P(|X-E(X)|>A)<^, 

where Y{Z) denotes the variance of Z. We start by computing the variance of 
X] E(X) has already been evaluated in (3). We have 



n 



= E + E inhh) - E(/.)E(/,)) (6) 
= nP(/i = 1)(1 - = 1)) + n{n - 1) {P(/i/2 = 1) - P'(/i = 1)} , 

so we focus on computing P(/i/2 = 1), i.e., the probability that both vertex 1 
and vertex 2 have odd stacks of pebbles. It is not too hard to verify that 



= ^ (Si - S2 - S3 + S4) , 



where 



and 



so that 



and 



^ \r,s,t ~ r - sj \nj \nj \ 
\r,s,t — r — sj \nj \ nj \ 
7^ yr, s, t — r — sy V nJ \nj \ n 

^ \r,s,t — r — sj \ nJ \ nJ \ n, 



Cov(/i,/2)=E(/i/2)-E(/i)E(/2) 
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It follows from (6) that 



We are now ready to state 

Theorem 2 Consider t distinct pebbles that are thrown onto the vertices of 
the complete graph Kn on n vertices according to the Maxwell Boltzmann dis- 
tribution. Set Aq = 1.5238 .... Then 



t — Aqh + ip{n)y/n =^ F{Kn is cover solvable) — >• 1 (n — >• oo) 



and 



t — Aqu — ip{n)y/n =^ E{Kn is cover solvable) ^0 (n — > oo), 

where (p{n) ^ oo is arbitrary. 

Proof Assume first that t — Aqu + (p{n)^. Then 



P(X >2n-t) 



= P 
= P 
> P 
= P 



> P 



X - E(X) >2n-t-'-\l - \ l- 



> 



n 
2 



X - E(X) > -n - Aon + ^ ( 1 



2x t 
n. 

n 



n 



n 



2 \ Aon+<^(n)-/n 



X-E(X)>--e-". + -^l--j 



- <^(n)Vnj 
-1 -'P{n)Vn ^ 



n 



X - E(X) > -e 



-2Ao 



exp 



^2ip{n) 



(p{n)y/n 



n 



X - E{X) > '-^e-^^° 



2^{n) 



;i + o(l))-(^(n)xA^ 



P (X - E(X) > -^{n)^l^ (l + e-'^°) (1 + o(l))) 
P (|X - E{X) I < (^(n) (l + e-'^«) (1 + o(l 



In (8), the first and second inequalities follow from the facts that \ —x< 
and e^^ — 1 < —x/ (1 + x) respectively. We next analyze the variance as given 
by (7). The first component ^(1 — (1 — ^)^*) of V(X) can easily be verified to 
be of the form |(l + iir(l + o(l))) (for some constant K) when t is as specified; 
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the second component 

n(n-l) if^ 4V r 2^] n(n - 1) [/ 4Y / 4^4^^^ 



4 |\ n J \ n J \ 4 |\ n / \ n 

on the other hand, may be bounded using the inequahties 

t{b - a)a*-^ < 6* - a* < t{b - a)h*-^ 
for 6 > a to yield, for t as above, 









nO ( 


V n 





n 



4 V? \ n \v? 

= m 

= e(n). 

It follows that N{X) = 0(n) and thus we have by (8) and Tchebychev's inequality, 
with K representing a generic constant, 

F(X >2n-t)>¥{\X - E(X)| < K ■ V^^(n)) > 1 - \ ^1 (n ^ oo), 
as asserted. 

The proof of the second half of the theorem is similar; we bound above instead of 
below to get, with t = Aqu — ip{n)y/n, 



F{X >2n-t) 

= P ( X - E(X) > + ^(^^^ ) 

+ ^ exp I (-2^0 + ^) (1 + o(l)) } + ^{n)^ 
= f[x- E(X) > ^e-2-4o jexp { + - l} + ^{n)^ 



<f[x- E(X) > |e-2^° . ^(1 + o(l)) + <^(n)v^) 
< P (|X - E(X)| > ^(n)v^ (1 + e-2^o(l + o(l)))) 

This completes the proof. □ 



10 



Remarks: Note that Y{X) = 0(n) implies that X is concentrated w.h.p. in an 
interval of length around K(X) = B(n). This is because of the inequality 

Y{X) < Range2(X)/4. 



4 Bose Einstein Cover Pebbling Threshold 



4-1 Exact Distributions 



In the Maxwell Boltzmann scheme, it is extremely difficult (though not impossible) 
to calculate F{X = x) exactly. Moreover, the formula for F{X = x) obtained thus is 
quite intractable. Surprisingly, this is not the case when we consider Bose Einstein 
statistics; we can derive such a formula by counting the number of configurations 
of size t on Kn with x odd stacks. Call this value (p{x,t,n). Clearly if x and t have 
different parity, or x > i or a; > n, (p{x, t, n) = 0. Suppose this is not the case. Then, 
we can find the configurations with x odd stacks of pebbles by placing one pebble 
on the X vertices that are to have odd stacks on them, and then distributing the 
remaining t — x pebbles on the n vertices of G in indistinguishable pairs. Thus, 
since the vertices with odd stacks may be chosen in ways, we have proved 

Proposition 3Ift and x have the same parity, and if x < mm{t,n}, then 



F{X 



In fact, the above exact distribution was used fruitfully in [25] to prove a weaker 
version of the main result of the next section. We choose, however, to prove a 
cover pebbling threshold for Bose Einstein pebbling by using more contemporary 
probabilistic tools. 



4-2 Poly a Sampling and Azuma's Inequality Yield Dividends 



There is a natural and sequential probabilistic process associated with Maxwell 
Boltzmann pebbling. We simply take t pebbles (balls) and throw them one by one 
onto (into) n vertices (urns) in the "natural" way that inspires many elementary 
problems in discrete probability texts. By contrast, the "global" Bose-Einsteinian 
positioning of t indistinct balls into n distinct urns - so that we obtain ("^^1^^) 
equiprobable configurations - does not appear to have a sequential process associ- 
ated with it. But it does. We first rephrase the problem - not as one associated 
with throwing balls into boxes but, conversely, as a sampling problem, i.e., drawing 
balls from boxes. In this light. Maxwell Boltzmann pebbling consists of drawing t 
balls "with replacement" from a box containing one ball of each of n colors, with 
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the understanding that the number of balls of color j drawn in the altered model 
equals the number of pebbles that are tossed onto vertex j d la the balls-in-boxes 
model. Bose Einstein pebbling can be recast in a similar fashion, but one needs to 
employ a process called Polya sampling. Polya sampling (or the Polya urn model) 
is described in [11] as a means of modelling contagious diseases and takes place as 
follows: Initially the urn contains one ball of each of n colors. After each draw, the 
selected ball is replaced together with another ball of the same color. In this mode 
of sampling, we lose the independence inherent to the with-replacement procedure, 
and, as a matter of fact, the selection process is not even Markovian but are able 
to "see" the sequentiality that will be critical in the sequel. As before, the number 
of times that color j is drawn can be set to equal the number of pebbles on vertex 
j, but do these two procedures yield the same probability model? We claim so, and 
here is a proof of this rather well-known fact: 

Lemma 4 Let Xj be the number of times the color j is drawn among the t draws. 
Then for any xi,X2, . . . ,Xn with ^ Xj = t, 

P(Xi =Xi,X2=X2,...,Xn = Xn)= fn+t-l\ 

\ t ) 



Proof First let us find the probability that the stated outcome appears in the 
following order: color 1 is first drawn xi times, then color 2 is drawn X2 times, etc., 
until we finally draw color n the last x„ times. Call this ordered event A. It is easily 
seen that 



p,^. ^ (1-2 xi)(l-2 X2)---(l-2 Xn) 

^ ^ n{n + l)---{n + t-l) 

_ XilX2l---Xnl 

~ n{n + 1) • • • (n + 1 - 1) ■ 

Moreover, a little reflection reveals that this probability is the same regardless of 
the order in which the balls are drawn. The total number of ways of ordering our 
configuration turns out to het\/x\\x2\ - ■ ■ Xn\- Thus, we can write the probability of 
obtaining a given configuration as 

f{Xi =Xi,X2=X2,...,Xn = Xn) = —. r • 

Xi\x2' ■ ■ ■ Xn' n\n + Ij • • • (n + r — 1) 
1 

This concludes the proof. □ 

We next provide the reader with background concerning the Azuma martingale 
inequality. Consider the following generic set up: {VL,T,f) is a probability space 
and (Yn) a sequence of random variables on (O, J^, P) that may not, in general, 
be independent. In our case, the YjS are the sequence of draws made in the Polya 
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urn scheme associated with the random cover pebbhng problem. Let X = Xt = 
X{Yi, . . . , Yf) be an objective function (in our case X is the number of odd stacks of 
pebbles), and consider the filtration (sequence of sigma algebras) {J^n), -^o = {Oj 
J^i = o"(Yi, . . . , Fj). Let EjX denote the conditional expectation of X with respect 
to Ti (with E = Eq) and set di = EjX — Ej_iX. Then it is well known that 
{di,J^i) is a martingale difference sequence, and that we have X-E(X) = 111=1 di- 
A key method used towards gaining an understanding of the concentration of X 
around K{X) is the method of bounded differences, also known as the Azuma (1967)- 
Hoeffding (1963) inequality [21]: 

Lemma 5 (Azuma-Hoejfding) For all A > 0, 

¥{\X - E{X)\ > A) < 2exp |-^^^^| . (9) 

where ||^||oo = (ess) sup|Z(a;)|. (9) may be made further applicable as follows: 
Letting {¥■*) be an independent copy of (Yi), we have 

E,_iX(Fi, ...,Yt)= EiX{Yi, Yi_i,Y*,Yi+i ...,Yt), 

so that di can be written as a single conditional expectation as follows: 

di = Ei{X{Yi, ...,Yt)- X{Yi, Y,_i,Y*,Yi+i ...,Yt)), 

and thus 



lldilloo = WMXiYi,. ..,Yt)- XiYi, Yi_i,Y*,Y,+i . . . ,yt)}||oo 

< ||X(Fi, ...,Yt)- XiYi, Y,^i,Y*,Yi+i . . . , Yt)\\^. (10) 

The philosophy behind the method of bounded differences is thus that the change 
in the value of X resulting from a change in a single input variable should be small. 
Moreover the bound in (10) shows us, after a moment's reflection, that we have, 
for our problem, ||(ii||oo < 2. Furthermore, Lemma 5 and (10) yield the following 
concentration for the number X of odd stacks in the Bose-Einstein scheme: 

¥{\X - E{X)\ > A) < 2exp |-^| , (Az) 

so that X is concentrated in an interval of length ^/nip{n) around E(X) whenever 
t ^ Kn. We are now ready to prove the main result of this section - one that 
features the golden ratio 7 = (1 + -\/5)/2: 

Theorem 6 Consider t indistinguishable pebbles that are placed on the vertices of 
the complete graph Kn according to the Bose Einstein distribution. Then, with 7 
representing the golden ratio (1 + •\/5)/2, 

t = ^n + ip{n)\/n F{Kn is cover solvable) — 1 (n — 00) 

and 

t = jn — ^p(ji)\fn P(-fCn is cover solvable) — (n — 00), 
where (p{n) — >■ 00 is arbitrary. 
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Proof We start by establishing tight bounds on E(X); unfortunately, the sum 

E{X) = nF{h = 1) 

does not appear to be tractable. Consider first an upper bound on P(7i = 1): We 
have, with a„ and representing, respectively, generic functions of the form 0{l/n) 
and 0(logra/n) whose exact form may vary from line to line. 



P(7i = 1) = P(vertex 1 has an odd number of pebbles) 

E \ n-2 ) 
(n+t-l\ 
j=l \ n-1 ) 

°° n-1 t t-1 t-2 t-2j + 2 



E n — J. 
. , n + t-V 



- ■ - -n + t-2n + t-^n + t-A n + t-2j 



< (1 + Q!n)> ■■■ 

- n^tn^t^ ' j^^n + t-3n + t-4 n + t-2j 

.7=1 ^ ' 



(1 + an) 



nt 



jn + t)^ ( t-1 V 

^ \n+t-3 J 



(1 + a„) 



= (i + «„);^. (11) 



Next note that for R to be determined, 



p(A = i)>f;. "-^ ' '-^ + ' 



^n + t- ln + t- 2n + t- 3n + t- 4 n + t-2j 
^^^^^' + ^-^2L^[n + t-2R 



2R 



1 _ / f-2R+2 y 
nt , \n+t-2Rl 

2 + 7 TT- (12) 



(n + ty , _ / t-2R+2 \ 

^ \n+t-2Rj 



We now pick i? in (12) so that 



t-2R + 2^'^^ 



'-[ n + t-2R 
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this may be done, e.g., with R = K log n, since in the range of f's that we are deahng 
with (namely i = 771 ± ip(n)^/n), we have 



if R is an appropriate multiple of logn. Equation (12) thus yields 



{n + ty^ "'(n-2)(n + 2t-4i? + 2) 
= (1 + «„)(! + 4);^ 

= (i + «^T2t- 

Equations (11) and (13) thus give 

and we are ready to derive our cover pebbling threshold. We first set t = ^n+ip{n)^/n 
to get, for some (and this is important) positive A, 



¥{X >2n-t) 

( nt 

>F[X- E(X) >2n-t (1 + . 

- V n + 2V 



X - E(X) > (2 - 7)n - ip{n 



)Vn . , ^. {l+lnj 

2 /' 1 I ip(n) 



X - E(X) > (2 - 7)n - ^(n) - ,(,'Mn) + 

n(l + 27) + 



> P^X - E{X) > (2 - 7)n - ^p{n)^/n 
Now (2 - 7)n - 7n/(l + 27) = 0, so (15) yields 



(15) 



F{X >2n-t) > P(X-E(X) > -ip{n)^ ^^v?(n)^/n + O(logn) ) 

V 1 + 27 J 

>P(X-E(X) > -B(^(n)Vn + 0(logn)) 

> P(|X-E(X)| < 593(n)Vn + 0(logn)) 

> 1 - 2exp [_&^^) + 0(logn))2 



8t 

0, ^ (16) 
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by Azuma's inequality as given by Equation (Az). 
Conversely, for t = ^n — ip{n)^/n, we have for some C > 0, 



P(X >2n-t) 

< P ( X - E(X) > (2 - 7)n + ^(n) 




n + 2(7n — 99(72) -^/n) 



n(7n — V9(n)y^) 





= P X - E(X) > (2 - 7)n + ^(n) 




(17) 



again by (Az). Equations (16) and (17) complete the proof. 



□ 



5 NP-Completeness of the Cover Pebbling Problem 

One of the obvious open problems that can be formulated as a result of our work is 
the following: What are cover pebbling thresholds for families of graphs other than 
It would certainly advance the theory of cover pebbling if one could uncover 
a host of results similar, e.g. to Theorems 2 and 6. Such results would provide a 
nice complement to those in [8]. Our results in this section show, however, that this 
task might not be as easy as one might imagine. Necessary and sufficient conditions 
for the cover solvability of a graph are likely to be complicated, and the best hope 
might thus be to establish necessary conditions and sufficient conditions that are 
not too far apart. 

The normal way to formalize the concept of the difficulty of a problem is to use the 
concept of compTitational complexity. Formally, we imagine a decision problem to be 
a set of infinite strings of characters (like data represented by bits in a computer.) 
A decision problem is said to accept a string if this set contains the string. Usually, 
we look for the best possible asymptotic upper bound (in terms of the length of 
the string) for the number of steps the fastest possible algorithm takes to deter- 
mine whether a given string is in the set. Informally, we think of decision problems 
being yes-no questions about a property of some class of finite mathematical struc- 
tures (graphs, integer matrices, etc.) and we ask how fast it is possible to correctly 
determine the yes or no answer in terms of the size of the input. 

For instance, some problems can be solved by an algorithm which takes only a 
number of steps which is bounded by a polynomial in the size of the input, while 
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others take at least an exponential amount of time to solve. The former class of 
decision problems is called P for "polynomial." The class NP, for "nondeterministic 
polynomial" is a bit more complicated; roughly speaking it is the set of decision 
problems for which a "yes" answer can be "checked" in polynomial time, given an 
appropriate piece of information. That is, if we call the class of inputs to the decision 
problem X, and the class of inputs which the decision problem accepts X', there 
exists a class Y of objects (called the certificates) and a function A : X xY ^ {0,1} 
which is computable in polynomial time, such that for any instance x E X oi 
the decision problem, there exists a y E Y such that A{x, y) = 1 if and only if 
X G X' . For instance, the decision problem which asks whether a given number is 
composite is easily seen to be in NP, because the composite numbers are exactly 
those with nontrivial divisor, and, given two numbers, it is easy to determine by 
division whether one is a divisor of the other. Also, any problem in P is also in 
NP, because any polynomial-time method of solving a problem is trivially also a 
polynomial-time method of verifying a yes answer. However, it is a celebrated open 
problem if the converse also holds and P = NP. 

Within NP, there is a class of problems, called the NP-complete problems, which is 
thought of being a set of problems which are at least as hard as any other problem 
in NP. This is because any instance a; of a problem D in NP can be translated by a 
polynomial-time algorithm to an instance x' of any A'P-complete problem D' such 
that x' is accepted by D' if and only if x is accepted by D. Therefore, if we could 
solve any A^P-complete problem in polynomial time, we could solve any problem 
in NP in polynomial time by translating it to an instance of this problem. Thus, 
the question of whether P = NP reduces to the question of whether any particular 
iVP-complete problem can be solved in polynomial time. 

We now show that the problem which asks if a configuration of pebbles on a graph 
is cover solvable is A^P-complete. It is worth noting that most complexity theorists 
speculate that P ^ NP, and therefore, when a problem is classified as A^P-complete, 
it is usually thought of as evidence of its difficulty. See [12] for a comprehensive 
theory of ATP-completeness, and Watson [23] and Milans et al. [19] for a more 
general exposition on the complexity of cover pebbling. 

Theorem 7 Let G be a graph, and C a confi.guration on G. Let \G\ = m and label 
the vertices of G as vi,V2, ■■ ■ Vm- Then C is cover solvable if and only if there exist 
integers n-ij > with I < i, j < m and riij = and riji = whenever {vi, Vj} ^ E{G) 
such that for all 1 < k < m, 

m m 

C{vk) + ^ njifc - 2 ^ nfe; > 1. 
1=1 1=1 

Proof First, suppose C is cover solvable. Then find some sequence of pebbling 
moves which cover solves C. Let Uij be the total number of pebbling moves from Vi 
to Vj in this sequence. Then after all the moves, there are exactly 

m m 

C{vk) + ^nik-2^ Uki 

1=1 1=1 
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pebbles left on Vk , which is always at least 1 because of the fact that this sequence 
of moves cover solves C. 

Conversely, suppose such numbers riij exist. This means that there does exist a 

sequence of moves that solves C, with riij moves being made from fj to Vj, but 
possibly with some illegal "negative pebbling" along the way. We show, however, 
that for each i,j it is possible to legally make riij moves from Vi to vj; since for each 
k, 

m m 

C{vk) + ^nik-2'Y^nki> 1, 
1=1 l=\ 

which leads to a cover solution of C. The main question thus is: In what order do 
we make these moves? We proceed in any arbitrary fashion, continuing to make 
pebbling moves as long as there exist vertices and Vji such that less than n^'j/ 
moves from to Vji have already been made and there are at least two pebbles 
on Vi'. If no such pair {vii,Vji} exists, then for each either nij moves have 

been made from i to j or else there is at most 1 pebble on vertex Vi. Let C be the 
configuration left on G after these moves and S be the set oi Vi & G for which the 
total number of moves from Vi is less than X]I=i 

If S" = then clearly for every 1 < i,j < m we have made n.^ moves from Vi to Vj 
and thus, for every k there are C{vk) + X]I=i ^ik ~ 2 Yll=i ^kl > 1 pebbles on Vk, so 
C'{v) > 1 for all V e G and we have cover solved C. 

If 5 7^ then consider the total number of moves that remain to be made from a 
vertex of S (the total X^jgg Z^JILi minus the number of moves that have already 
been made from vertices in S.) By the definition of S, this total is at least l^l, since 
at least one move remains to be made from every vertex in S. The total is exactly 
\S\ only if exactly one move remains to be made from each vertex. Also, C'{v) < 1 
for all V & S, for a total of at most IS'I pebbles. Consider the remaining moves, each 
of which must originate from S. Each of these moves, if executed, would remove one 
pebble from S if they also end at a vertex in S, and two if they end at a vertex 
outside of S. Thus each move must both begin and end in S, for otherwise S would 
be left with a negative total number of pebbles at the end of the pebbling sequence, 
which is impossible. Even if all moves begin and end in S, however, we end up with 
at most pebbles on S" - which too is impossible since we have assumed that the 
moves cover-solve, and so there must be one pebble on each vertex of S. Thus we 
must have ^ = and we have cover solved C. □ 

Corollary 8 The cover solvability decision problem which accepts pairs {G, C} if 
and only if G is a graph and C is a configuration which is cover solvable on G is in 
NP. 

Proof The above theorem gives the appropriate certificate of cover solvability, any 
list of integers Uij which satisfy the equation in Theorem 7. Indeed, Theorem 7 shows 
that cover pebbling is equivalent to a special case of the A'P-complctc problem of 
integer programming, which asks, given an n x m integer matrix A and an n- 
dimensional integer vector h if there exists an m-dimensional integer vector x such 
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that Ax > b, holds componentwise. Having reduced cover solvabihty to a special 
case of this NP problem, we know that cover solvability is also in NP. 

We now pause to point out another corollary which will be needed later and which 
is interesting in its own right: 

Corollary 9 Let G be a graph, C a configuration on G. If the sequence of pebbling 
moves Q = {qi,q2 ■ ■ ■ Qk) cover solves C, and it is possible to make the sequence 
Q' = {qii , qi2 , ■ ■ ■ qii) of moves (with 1 < ij < k for all j but with no particular 
requirement on the order of the ij,) then the configuration C' obtained from C after 
the moves {qii, ! ■ ■ ■ Qii) cover solvable. 

Proof The point of this corollary is that the order of our pebbling moves can't 
matter. To show this, we simply note that if it were possible to somehow execute 
the remaining moves from Q which are not in Q' , they would solve C'. By Theorem 
7, it is thus possible to solve G' . □ 

Now we turn our attention to showing that the cover solvability decision problem 
is iVP-hard, that is, that any instance of any problem in NP can be translated to 
an instance of cover solvability in polynomial time. The usual method of showing 
that a problem A is A'P-hard is to find an A/P-complete problem B for which any 
instance of B can be translated into an instance of A in polynomial time. Then for 
any instance of any problem in NP we can translate it in polynomial time to an 
instance of B, then translate this instance into an instance of A. For cover solvability, 
we will use a known AfP-complete problem known as "exact cover by 4-sets." Indeed, 
the corresponding and seemingly simpler problem of perfect cover by 3-sets is also 
iVP-complete, but for our purposes, the 4-set problem is more useful. 

Theorem 10 (Karp) [17] Let the exact cover by 4-sets problem be the decision 
problem which takes as input a set S with 4n elements and a class A of at least n 

^-element subsets of S, accepting such a pair if there exists an A' A such that A' 
is a class of disjoint subsets which make a partition of S, that is they are n subsets 
containing every element of S. This problem is NP-complete. 

Theorem 11 The cover solvability decision problem is NP-complete 

Proof Having shown that this decision problem is in NP, it remains to be shown 
that it is NP-haid. Given a set 5 = {si, S2, ■ ■ ■ s^n} and a class A = {ai, 02, . . . a^} 
of four element subsets of S, that is, an instance of the exact cover by 4-sets problem, 
we construct a graph G' and a configuration G' on G' in the following manner: We 
create a set of vertices T = {ti,t2, ■ ■ ■ t4n} which will be thought of as corresponding 
to the elements of S, and a set of vertices B = {61, 62, • • • &m} which will be thought 
of as corresponding to the members of A. Let G'{t) = for all t G T and let G'{b) = 9 
for all b £ B. We make edges between B and T in the intuitive way, including {bi, tj} 
if tj £ bi . Additionally, create a vertex v and a path of length m — n which has one 
terminal vertex v and the other called w. Let G'{v) = 2™~" — (m — n) + 1, G'{w) = 
and G'{u) = 1 for all u between v and w on the path. Finally, create vertex classes 



19 



B' = {b[,b'^...b'„,} and B" = {b'i,h'i . . X.}. creating edges {6,, 6^}, {b',X} and 
{b'l, v] for all i. Let C'{u) = 1 for all « € 5' U 5". (Figure 1.) 



ti t2 ts ti ts ^7 ^8 




w9 

Fig. 1. A graph that corresponds to the exact cover by four 4-sets problem, 

ai = {Sl, S2,S3, S4}, 02 = {S3, S4, S5, Sq}, 03 = {55, Sq, S7, Ss}. 

Clearly, this construction can be made in polynomial time in the size of the pair 
{S,A}. Indeed, we have made only 5\A\ + 1 vertices and 8\A\ — \S\ edges. In order 
to finish the proof, we claim that C is solvable if and only if A contains a perfect 
cover of S. 

First STipposc that A contains a perfect cover A' = {ajj, ai,^, . . . a.j^} of S. Then for 
each vertex in B which is a bi- for some 1 < j < n, we use 8 of the pebbles on 
this bi-, two each to cover the four vertices of T to which it is adjacent. Because 
of the fact that A' is a perfect cover and the way we constructed G' , we now have 
exactly one pebble on every vertex of T. Furthermore, we have m — n vertices in 
B which still have 9 pebbles each on them. Because v is at distance 3 from each of 
these vertices, we can use 8 pebbles from each of these vertices to move one pebble 
each onto v from these m — n vertices. This leaves 2"^""' + 1 pebbles on u, which 
is enough to move one pebble onto w while leaving one pebble on v. After this is 
done, we have exactly one pebble on every vertex of G' , and we thereby know that 
C is solvable. 

To show the converse, suppose that A does not contain a perfect cover of S. Suppose 
as well that C is solvable on G' . Clearly, the sequence of pebbling moves which solves 
C must contain (at least) one move to t for every t eT. Clearly, each of these moves 
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must originate from B, and no more than 4 can originate from any one vertex of 
B. Since A does not contain a perfect cover of S, it cannot be the case that these 
moves originate from exactly n vertices in B. 

We make these 4n moves immediately from C, using Corollary 9 to sec that the 
resulting configuration must be solvable (we use the fact that no more than 4 of 
these moves can originate from any one vertex in B to see that it is indeed possible 
to make these moves). In addition to the one pebble left on every vertex of B to 
ensure they remain covered, there are now 8(m — n) pebbles on B, but they are 
not in m — n groups of 8 pebbles because the moves we made originated from more 
than n vertices oi B. In order to reach w, we clearly need to move m — n pebbles 
onto V while leaving the rest of the graph covered. Clearly, this is only possible if 
all 8(m — n) extra pebbles are moved by a path of length 3 onto v. However, only 
one such path is available for any group of these pebbles. But any group of less 
than 8 pebbles cannot increase the number of pebbles on v by moving along this 
path, while leaving the vertices of the path covered. Since there are not indeed m — n 
groups of 8 pebbles on B we see that it is impossible to gather the pebbles necessary 
to reach w and thus our configuration is not solvable, which is a contradiction. 
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